function read_marc,file

a=read_ascii(file)
b=reform(a.field1(0,*))
c=reform(a.field1(1,*))
d=reform(a.field1(2,*))

h={x:b,y:d,e:c}

return,h
end

pro plt1,h

nn=size(h,/n_elements)
lab=['Z 0.01-0.10','Z 0.10-10','Z 0.01-10.0','','','']
for i=0,nn-1 do begin
	k=read_marc(h(i))
	if( i eq 0) then plot,k.x,k.y,/xlog,/nodata,yrange=[0.5,1.5],ytitle='MARC CORR',Xtitle='R [Mpc/h]'
	oplot,k.x,k.y,linestyle=i,color=mycolor(i),psym=5
	errplot,k.x,k.y-k.e,k.y+k.e,linestyle=i,color=mycolor(i)
	xyouts,10.02,i*[1.5-1.0]/nn+1.05,lab(i),col=mycolor(i)
endfor
oplot,[0.001,1000.0],[1.0,1.0],linestyle=8

end

pro plterr,x_tab,u_tab,u_stab,col,gg
ntab=size(x_tab)
Ntab=ntab(1)
dxe=float(gg)/float(ntab)/2
FOR i=0,ntab-1 DO BEGIN
      oplot,[x_tab(i),x_tab(i)],u_tab(i)+u_stab(i)*[-0.5,0.5],color=mycolor(col)
      oplot,x_tab(i)+dxe*[-1,1],u_tab(i)+u_stab(i)*[0.5,0.5],color=mycolor(col)
      oplot,x_tab(i)+dxe*[-1,1],u_tab(i)-u_stab(i)*[0.5,0.5],color=mycolor(col)
END
end
